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Abstract 

In nonspherical nuclei at J = J c +0 the relationship between the angular momentum 
and angular velocity immediately after backbending is the same as in the limiting case 
J — J c — > oo. This indicates that there is a unique type of cancellation of the deviations 
from a rigid-body moment of inertia in the upper phase J > J c . An integral relationship 
is found which expresses this cancellation quantitatively. This formula permits J c to 
be calculated for the rotational bands of the even-even nuclei studied and the results 
are in agreement with those obtained by other methods of locating the Curie point. 
For the ground state band of W 170 the cancellation of the reciprocals of the true and 
rigid-body moments of inertia can be verified directly. The condition for the stability of 
the rotation of a nonspherical nucleus is analyzed in the Appendix in close connection 
with the problem of a reasonable definition of the concept of a variable moment of 
inertia. 



1 Introduction 



Even before the discovery in 1971 of the singularity exhibited by the rotational band at 
a certain critical value of the nuclear spin J = J c , it was gradually becoming clear that 
we are actually dealing with a situation where the moment of inertia J of a nonspherical 
nucleus varies most significantly. In choosing a reasonable definition of this concept it is 
desirable to keep in mind the following considerations, in addition to purely aesthetic ones. 
First of all, for completeness and internal consistency of the theory it is important that both 
definitions of the moment of inertia, that via the Lagrangian and that via the Hamiltonian, 
be equivalent. In addition, in accordance with general physical considerations it is natural 
to expect that precisely at negative values of the correctly defined moment of inertia the 
rotation becomes unstable; this is discussed in more detail in the Appendix. 

*Zh. Eksp. Teor. Fiz. 76, 1506-1514 (1979) [Sov. Phys. JETP 49, No. 5, 765-769 (1979)] 
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The following definition meets all these requirements: 

h0 _dE h 2 _d{htt) _d 2 E 

(the notation is the same as in Ref. 1). We emphasize that these prerequisites do not create 
any practical difficulties in a concrete comparison with experiment since in formulas like 
(1) it is possible to replace the derivatives with respect to J by the corresponding directly 
observable finite differences. In fact, the widely known Bohr-Mottelson formula 

E^^J(J + 1) (2) 

has the property that the finite differences calculated with it agree with the result of the 
formal differentiation according to (f). On the other hand, for J ^> 1 the replacement of the 
derivatives by finite differences suggests itself automatically and no special problems appear. 

In our earlier study [1] this phenomenon was viewed as a smooth, continuous rearrange- 
ment of the angular momentum coupling scheme in the nucleus. For adiabatically slow 
rotation, J <C J c , the internal state of the system is formed mainly by a "nucleon-nuclear 
symmetry axis" type of interaction, which is due to the axially symmetric deformation of a 
nucleus "at rest". However, in reality the nucleus is rotating, so that the nucleon-rotation 
axis interaction (which here emerges as a certain non-adiabatic correction) is always effective 
to some degree. In the entire region J < J c there is a complicated "nucleon-nuclear sym- 
metry axis" plus "nucleon-rotation axis" coupling scheme of relatively low symmetry and 
the two interactions, generally speaking, are comparable. However, in a sufficiently strong 
rotation "field" the mechanical angular momenta of the individual quasiparticles are aligned 
parallel to the fi|| J direction and cease to be oriented along the vector n. It can be said that 
this corresponds to the simplest, most symmetric nucleon-rotation axis coupling scheme, 
not directly affecting the direction n of the nuclear axis. Then the vector n remains "free," 
that is, it is actually distributed isotropically for J > J c . As a consequence of the increased 
symmetry of the rotational state at the point J = J c the angular velocity of the rotation falls 
abruptly by some amount A(hQ) and in the isotropic upper phase the moment of inertia 
displays the seemingly paradoxical limiting behavior 

I = j/(J-J c ), J-J c «i// (3) 

(j is some coefficient depending on the particular nucleus and I is the rigid body moment 
of inertia). 

The theory that we developed earlier [1] does not claim to be able to calculate the specific 
values of such parameters as J c , the discontinuity A(hQ), or the coefficient j for individual 
nuclei. In relation to this it should also be noted that the less symmetric lower phase is 
considerably more difficult to study theoretically. The main result here is the square-root 
law 

Q (x (J c - J) 1/2 , (4) 

according to which the static quadrupole moment vanishes near the Curie point. It seems 
almost certain that it is not possible to quantitatively determine the value of the critical 
angular velocity fi mc , for example, from only deductive considerations. However, it is re- 
markable that for a given J c the value of Q nc to which the rotational velocity falls after its 
discontinuous decrease, can be calculated in a closed form. 
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2 The angular momentum and angular velocity in the 
upper phase J > J c 

After reaching the rigid body value of the moment of inertia 

/ = J , J - J c > j/I , (5) 

the rotational motion of the two components of nuclear matter can be assumed to be fully 
concurrent. Let us analyze this rotation in detail, appealing to a physical manifestation of 
the proton component like the magnetic moment. It is obvious that here 

9 = 9o = j (6) 

is valid for the gyromagnetic factor. 1 

However, the nuclear magnetic moment can be calculated in a somewhat different way, 
using the Larmor theorem (see Ref. [3], for example) in, so to speak, its differential form. 
If it is favorable for a spherical nucleus to rotate in the above manner, then in an applied 
magnetic field H it can be viewed as rotating, but with an angular velocity decreased by the 
Larmor value 

n* = (7) 

2m p c 

(m p is the proton mass). Therefore, the rotational properties of the system that we are inter- 
ested in here are described by the function E(Q), while the change of energy is —Q L dE/dQ. 
According to the usual point of view, the nucleus at a level J is a particle having magnetic 
moment \i and additional energy —fiH in the field (both the angular momentum vector and 
the magnetic field vector are assumed to be directed along the z axis). 

Let us equate the two expressions for the additional energy of the system in the magnetic 
field: 

-%n L = (s, 

Furthermore, transforming the derivatives according to (1) and substituting formula (7), we 
also take into account the fact that according to the conventional definition of the nuclear g 
factor, the nuclear magneton eh/2m p c must serve as the unit of measurement of the magnetic 
moment: 

g = 9o m/hJ. (9) 
Finally, substituting here the values (5) and (6), we find 

h.j = i n, j-j c >j/7 . (io) 

This seemingly natural recovery of the simple proportionality between the angular mo- 
mentum and the angular velocity sheds additional light on the physical nature of the upper 

1 lt is well known thai in adiabatic region the nuclear matter is taken into rotation only partially and 
the proton component is involved into rotation less than the neutron component. For the lowest levels of 
even-even nuclei it is observed experimentally that g < go; see, e.g. [2]. This physical non-equivalence of 
protons and neutrons must disappear in the solid state limit of rotation. It is not necessary to take into 
account spin magnetism of separate nucleons in our macroscopic approach. 
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phase. For J — J c ^> j / 1 its characteristics cannot depend on such non-universal parameters 
as the coefficient j, for example. Roughly speaking, here we have a fairly clear interpretation 
of the ordinary rotation of a rigid body. However, if we move down the band in the negative 
direction of the J axis the stability of this regime deteriorates. In close relation to this the 
moment of inertia undergoes an odd, non-monotonic change. First it falls off, then passes 
through a minimum, and then approaches the pole according to (3). In the region near the 
transition point J — J c ^ j / Iq of the variable moment of inertia the rotation, of course, can 
no longer be viewed as simply rigid body rotation. 

Let us return for a while to the region where the moment of inertia is constant I = Jo- 
in our earlier study [1] in a somewhat formal manner we found 

hVt = hVt nc + h 2 (J -J c )/I (11) 

for the rotational velocity. Now, comparing this to formula (10), we have 

hj c = i n nc . (12) 

Therefore, in reaching the pole J = J c + the effect of the above-mentioned non- 
monotonic change of the moment of inertia on the rotational velocity is cancelled. It is 
easy to find an integral relation expressing this cancellation in a quantitative form. For 
this according to the second formula in (1) we integrate the inverse value of the moment of 
inertia from J c up to some large J; then, performing an even more trivial integration of the 
analogous expression with the rigid body moment of inertia Iq, we take into account relation 
(12) and the fact that for J — > oo equation (10) is valid. As a result we easily find 

---)dJ = 0. (13) 

For practical purposes it is convenient to write this relation in a dimensionless form: 

PCX) 

/ (l-I /I)dJ = 0. (14) 

JJc 

In addition to the existing means of determining the location of the phase transition 
point, formula (12) permits J c to be calculated according to the observed velocities Q nc . In 
Table I we compare the results of processing the experimental data by different methods 
for the twenty-eight rotational bands of nonspherical even-even nuclei that we studied. The 
/^-vibrational bands of gadolinium and dysprosium are starred. In column 2 we indicate, of 
necessity tentatively, the location of the discontinuous decrease of the angular velocity of 
rotation. The value J^. xtr is the result of extrapolation according to formula (3) (see Ref. 1 
for more details). The last column gives the value of the critical spin calculated according 
to (12). In calculating the rigid body moment of inertia 1$ we used the value of the nuclear 
radius obtained in Ref. 1 r = 1.1 x 10~ 13 cm, which is also in agreement with the data on 
electron scattering. It corresponds to the working formula 

ft 2 85900 n wl 



Jc 
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Table I 



Nucleus 
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jextr 


TTT 

jtneor 


Nucleus 


jexp 


jextr 
J c 


TXT 

jtneor 


56Bagg 4 


> 11 




< 13.4 


771158 

68^90 


13 - 


15 




12.6 


"Ro 126 

56-Da 70 


11 - 13 




12.5 




25 - 


27 


25.9 


23.1 


p„128 
58^ e 70 


11 - 13 




10.6 


"I7V160 

68-CJ92 


13 - 


15 




14.7 


58 Ce^ 


9-11 


11.0 


9.7 


68- t!jr 94 


13 - 


15 


15.7 


15.2 


58 Ce^ 2 


11 - 13 




10.3 


68^ r 96 


15 - 


17 




14.4 


58^76 


9-11 




9.5 


"\/"U^ 164 

70 1 b 94 


13 - 


15 


13.0 


14.0 


64 Gd^ 4 


> 17 




< 15.9 


70 Ybg^ 6 


13 - 


15 


15.5 


14.5 


(~~\ A 154* 

64^%) 


11 - 13 


12.4 


10.5 


"VU168 
70 I D98 


15 - 


17 




18.3 


66Dyg| 4 


13- 15 




14.1 


70 I D 100 


15 - 


17 




18.6 


156 

66JJy 90 


15 - 17 




16.0 


TTX168 

72nt 96 


13 - 


15 




13.5 


66-^90 


11 - 13 




9.5 


72Hf 9 |° 


15 - 


17 




18.6 


66^92 


15 


14.8 


16.6 


-117-170 
74 W 96 


11 - 


13 


13.6 


13.1 


D„160 
66^94 


15-17 




16.0 


76^S 106 


13 - 


15 


11.9 


16.3 


T7V156 

68^88 


11 - 13 


13.2 


13.8 


76^Sl08 


13 - 


15 


15.1 


18.4 










76^S no 


> 15 




< 17.4 



On the whole the agreement appears to be satisfactory. It should, however, be noted that 
the nearness of the magic number N = 82 can sometimes manifest itself in very unexpected 
ways. It can be supposed that the second singularity found experimentally in the 
ground state rotational band of Er 158 (Ref. 4) is due to this; see also Table I. In our opinion 
the construction of more detailed hypotheses on the nature of the "intermediate phase" 
Jc 1 ^ < J < Jc 2 ^ is still somewhat premature. Intuition suggests that as the magic nucleus is 
approached the probability of similar surprises increases. In fact, for Er 156 we already have 
circumstantial evidence of this supposition. According to the data of Ref. 5, after the first 
phase transition the moment of inertia, as usual, passes through a minimum, but then rises 
very steeply, reaching / = 1.38J for J = 22. The available information indicates that Er 156 
will also have a second backbending region. 

As an illustration of the agreement between the different methods of finding J c it is worth 
noting that the rotational velocity Kl nc varies within a fairly wide range: from 216 keV for 
W 170 to 340 keV in the ground state rotational band of Ba 126 . The rigid body moment of 
inertia Iq increases by roughly a factor of two throughout the entire table. 

For nonspherical nuclei sufficiently far from the possible effect of the magic numbers 
(see above), experimental data on the upper phase are far from abundant. Therefore, the 
possibilities of directly verifying the integral relation (13) or (14) at the present time are 
very limited. Only for W 170 does the moment of inertia, after passing through the minimum, 
approach the rigid body value with an accuracy of about 4%. In Fig. 1 we give the graph 
of the function 1 — I /I in this case. The accuracy with which the integral of this function 
becomes zero can be considered satisfactory. 
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3 Conclusions 

The reason that for tungsten there are few experimental points on a large part of the moment 
of inertia curve in the upper phase is the following: the ratio j/I is not large, amounting only 
to 1.88. However, there are nuclei for which the value of this parameter is much larger. For 
Dy 158 , for example, the number of neutrons N = 92 is sufficiently far from the magic number 
and j / Iq = 6.76. Here further study of the upper phase makes it possible to construct a 
more accurate moment of inertia curve using a considerably larger number of experimental 
points. Beginning, for a rough estimate, from the assumption of similarity to the curve in 
Fig. 1, we conclude that in the ground state band of Dy 158 it is of interest to measure the 
location of the rotational levels up to J ~ 40, beyond which the moment of inertia becomes 
practically the rigid body value. 

In addition to the data of Ref. 6 cited in Ref. 1 on the radiation of Ce 134 , 2 verification 
of the "area theorem" (13) for different nuclei, including cases with j/I 3> 1, would give 
interesting new material for judging the validity of the theory treating backbending as a 
macroscopic quantum phenomenon. 

We are grateful to I.M. Pavlichenkov for discussing the results of this study. 



Appendix. The condition for rotational stability of a 
nonspherical nucleus 

The considerations discussed below apply equally to either phase J ^ J c . Let us view the 
level Jo of the ground state rotational band as minimizing the total energy of the nucleus 
for a given value of the conserved angular momentum of the entire system: 

£ = min, J = J . (A.l) 

Using the Lagrange multiplier method, we shall drop the auxiliary condition and uncondi- 
tionally require that 

E-XJ = min. (A.2) 

This can be rewritten as 

5(E - XJ) > 0. (A.3) 

Moving now from the minimum along the actually realized rotational band E(J), let us 
calculate the energy change SE with accuracy to second-order terms inclusive: 



The value 

A = (dE/dJ) J=Jc (A.5) 

of the Lagrange multiplier ensures the correct location of the extremum and is a minimum 
for 

< A6) 



2 Similar data were obtained in [7] for Ce 132 and Ce 130 . 
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Here we have taken into account the fact that the point J = J was chosen arbitrarily. 
Comparison with formula (1) gives 

/ > 0. (A.7) 

Therefore, the requirement that the moment of inertia be positive emerges here as the 
stability condition. 

The question of the rotational stability of a nonspherical nucleus can also be approached 
from a somewhat different viewpoint. Let us as usual denote the spherical angles giving the 
orientation of the vector n (if convenient it plays the role of the rotational variable, but no 
separate "rotational Hamiltonian" corresponds to it) in the stationary space by 9 and 0. 
The state of motion in the angle 9 is, in general, not "pure." The situation with the motion 
in the azimuthal angle is different: because of the conservation of the z component of the 
total angular momentum it corresponds to the separate wave function 

Vv ot = (2TxT x l\ m \ M = J z . (A.8) 

For J z = J ^> 1 the angular momentum is directed along the z axis and the situation is 
semiclassical: the fully described state (A.8) of regular precession can in this limit be viewed 
as changing into classical motion in a cyclic trajectory, corresponding to the variation of the 
azimuthal angle 0. 

In order to remove possible doubts about the validity of using a purely mechanical ap- 
proach here, we recall that the free rotation of a body in thermodynamic equilibrium is 
not accompanied by friction (see Ref. 8, for example). We shall proceed directly from the 
principle of least action 

f Ldt = min, 0(tx) = 0!, 0(t 2 ) = 2 , (A.9) 

(see Ref. 9, for example). The extreme values of the angles 0i and 2 are viewed as constants 
which are not varied: 

<ty(ti) = 8<t>{t 2 ) = 0. 

However, let us first consider how the moment of inertia is expressed in terms of the La- 
grangian L. We have 3 

dE 

E — E(M), j^ = ^ L = L(Q). 

Let us now express the Hamiltonian E(M) in terms of the Lagrangian: 

E = VLdL/dVL - L. (A.10) 

Differentiating (A. 10) once with respect to M and cancelling f2 from both sides, after some 
elementary manipulations we easily find 

d 2 L d 2 E _ 

dQ 2 dM 2 ~ ' 

3 Because of purely classical formulation of the problem we take M = hJ and use from now on usual units. 



7 



Comparing this to formula (1), we finally find 

/ = d 2 Ljd£l 2 = {d 2 E/dM 2 )-\ (A.f f) 

This is precisely the relation between the moment of inertia and the Lagrangian that we 
would naturally expect; see also the preliminary discussion at the beginning of this article. 
Let us now return to the principle (A. 9), rewriting it in the form 



5 f 2 L(Q)dt > 0. (A.12) 
Ju 



Free rotation can occur only as uniform rotation, that is, when the angular velocity is 
constant: 

n = (A.i3) 

The deviation 5<p(t) from this simple law of motion gives rise to a time-dependent variation 
of the velocity 

SQ(t) = d5(f)(t)/dt. 

Let us find the corresponding variation of the action with accuracy to second-order terms 
inclusive: 

It is obvious that 

f 5Q{t)dt = 5(j){t 2 ) - Scfrih) = 0, [ [5Q(t)] 2 dt > 0. 

Therefore 

5? > 0- (A.U) 

Comparing this with formula (A. 11), we find the stability condition (A. 7). States with a 
negative moment of inertia as violating the principle of least action are unstable and cannot 
be realized as stationary rotational energy levels of a nonspherical nucleus. 

In conclusion let us clarify the nature of the limiting behavior of the Lagrangian of the 
upper phase. It can be found explicitly in the region near the transition point, where the 
moment of inertia has a pole behavior (3). Substituting this formula into (A. 11), using formula 
(30) from Ref. 1 we can express the angular momentum in terms of the angular velocity and 
then integrate twice with respect to Q: 

L n (Q) = ^(2%') V2 (fi - Q nc f' 2 + I n nc n - E c . (A.15) 

The constants of integration are here chosen such that the first derivative gives the angular 
momentum M = hJ and the energy (A. 10) coincides for Q = Q nc with the nuclear excitation 
energy E c at the Curie point. We see that the Lagrangian of the upper phase has a singularity 
(a branch point) at the point of the phase transition Q = Q nc . 
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Figure caption 

Fig. 1. The integral relation (14) for the W 170 nucleus. The area under the "anomalous" 
(/ > Iq) part of the curve is cancelled with an accuracy of about 10%. 
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This figure "fig.gif" is available in "gif" format from: 



http://arXiv.org/ps/nucl-th/0405001vl 



